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Abstract 

Gamma ray burst (GRB) fireballs provide one of very few astrophysical environments where one can con- 
template the acceleration of cosmic rays to energies that exceed 10^° eV. The assumption that GRBs are 
the sources of the observed cosmic rays generates a calculable flux of neutrinos produced when the protons 
interact with fireball photons. With data taken during construction IceCube has already reached a sensitiv- 
ity to observe neutrinos produced in temporal coincidence with individual GRBs provided that they are the 
sources of the observed extra-galactic cosmic rays. We here point out that the GRB origin of cosmic rays is 
also challenged by the IceCube upper limit on a possible diffuse flux of cosmic neutrinos which should not be 
exceeded by the flux produced by all GRB over Hubble time. Our alternative approach has the advantage 
of directly relating the diffuse flux produced by all GRBs to measurements of the cosmic ray flux. It also 
generates both the neutrino flux produced by the sources and the associated cosmogenic neutrino flux in a 
synergetic way. 
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1. Motivation 

There is compelling evidence that the collapse of a massive star to a black hole is the primary engine 
of long-duration Gamma Ray Bursts (GRBs). The phenomenology that successfully accommodates the 
astronomical observations is that of the creation of a hot fireball of electrons, photons and protons that 
is initially opaque to radiation. The hot plasma therefore expands by radiation pressure and particles are 
accelerated to a Lorentz factor F that grows until the plasma becomes optically thin and produces the GRB 
display. From this point the fireball is coasting with a Lorentz factor that is constant and depends on its 
baryonic load. The baryonic component carries the bulk of the fireball's kinetic energy. The energetics 
and rapid time structure of the burst can be successfully explained by shocks generated in the expanding 
fireball [iHS]. Here, the temporal variation of the 7-ray burst of the order of milliseconds can be interpreted as 
the collision of internal shocks with a varying baryonic load leading to differences in the bulk Lorentz factor. 
Electrons accelerated by first order Fermi acceleration radiate synchrotron 7-rays in the strong internal 
magnetic field and thus produce spikes in the burst spectra of the order of seconds [H [S] . The collision of the 
fireball with interstellar gas forms external shocks that can explain the GRB afterglow ranging from X-ray 
to the optical [6l[7]. (See also [8] for a recent review on theoretical models.). 
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It has been pointed out that firebah baryons may be the source of uUra-high energy (UHE) cosmic rays 
(CRs) with energies extending to at least 3 x 10^" eV. Baryons are inevitably accelerated along with the 
electrons in the expanding fireball with their energies boosted by the bulk Lorentz factor. It has been shown 
that a typical GRB environment can naturally satisfy the requirements to produce UHE CRs, most likely 
protons. It is not easy to conceive of a mechanism where nuclei survive acceleration in a GRB fireball 
environment. (Heavy nuclei could be synthesized in magnetically-dominated jets as in the proto-magnetar 
model of GRBs, see e.g. [S]-) The results of CR observatories disagree with HiRes [TU] claiming protons and 
Auger heavier primaries. These conflicting results may just illustrate that, given the poor knowledge of 
hadronic interactions more than one order of magnitude above LHC energies, is not sufficiently known to 
derive a definite result |12j . 

Though GRBs satisfy the necessary conditions for accelerating protons to UHE, it is problematic how 
these protons may eventually be ejected as CRs: protons are magnetically confined to the expanding fireball 
and its adiabatic cooling will reduce the maximum proton energy significantly [13l[T4]. However, this does 
not concern neutrons that are frequently produced in p7-interactions of accelerated protons with fireball 
photons by processes like p7 — > A+ — > wk^ . Cosmic ray protons could thus be identified as neutrons 
from p7-interactions that can escape from the magnetic environment and /3-decay back to protons at a 
safe distance. A smoking-gun test of this scenario is the production of PeV neutrinos from the decay of 
the charged pions inevitably produced along with the neutrons |15H18j . This neutrino flux is in reach of 
large-scale neutrino telescopes like IceCube [El [20] . (Possible emission of neutrinos associated to the GRB 
progenitor and prior to the 7-ray burst have been discussed in |21) . Neutrinos from proton interactions in 
late internal or external shocks during the afterglow phase have been considered in [22n24) .) 



In this paper we calculate the diffuse flux of neutrinos produced in association with GRB cosmic rays 
by directly fitting the proton spectra from the decay of fireball neutrons to CR data from HiRes. As a 
byproduct we also obtain the flux of so-called GZK neutrinos produced when the cosmic rays interact with 
microwave and infro-red/optical background photons |251I26| . The accompanying photon flux resulting from 
the electromagnetic cascading of the neutral pions peaks in the MeV-GeV energy region; we verify that it 
does not exceed the extra-galactic diffuse 7-ray flux inferred by Fermi-LAT 27J. Our main conclusion is 
that the predicted flux exceeds the upper bound on a diffuse flux of cosmic neutrinos obtained by IceCube 
from a year of data taken with half the instrument during construction. Facing this negative conclusion, 
we subsequently investigate the dependence of the predicted neutrino flux on the cosmological evolution 
of the sources as well as on the parameters describing the fireball. Although the latter are constrained by 
the electromagnetic observation as well as by the the requirement that the fireball must accommodate the 
observed cosmic ray spectrum, the predictions can be stretched to the point that it will take 3 years of data 
with the now completed instrument to conclusively rule out the GRB origin of UHE CRs. 

We work throughout in natural Heaviside-Lorentz units with h ^ c ^ e^) ^ fi^ ~ 1, a — e^/(47r) ~ 1/137 
and 1 G ~ 1.95 x IQ-^eV^ 

2. Fireball Model 

In the fireball model the GRB central engine produces a heated optically thick plasma of leptons, photons 
and baryons which is initially at rest. The fireball expands adiabatically by radiation pressure until it 
becomes optically thin. From this point on the fireball is coasting with a Lorentz factor F which depends on 
its baryonic load. The expanding plasma flow is shocked producing shells with varying baryonic load; this 
results in internal shells with different velocities that in their collision produce internal shocks. Although in 
the baryon-rich flreballs that would produce EHE cosmic rays most of the kinetic energy of the flreball is 
carried by the baryons, the internal shocks will convert a fraction Ee — U'e/^icm to leptons and Eb = U'^/Ul--^j^ 
to magnetic flelds supported by the plasma. Here and in the following, primed quantities refer to values in 
the comoving plasma frame, whereas unprimed quantities are reserved for the observer's frame. 

Study of the emission of GRB has resulted in "benchmark" parameters describing GRB fireball with 
the efficiency for lepton acceleration from bulk energy taken to be 10%, i.e. Se — 0.1. The fraction of 
magnetic field energy and hence the value of the turbulent magnetic field strength can be inferred from the 
observed peak of the 7-ray synchrotron spectrum with Eb — 0.1, corresponding to equipartition of energy 
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in leptons and magnetic field. The "benchmark" parameters are required to produce the observed peak 
emission at MeV energies resulting from the synchrotron radiation of GRB electrons and must be consistent 
with afterglow observations. The typical radius of on internal shock is given by the speed-of-light distance 
implied by the duration of the spikes observed in the GRB display. Because the shells are boosted toward the 
observer by a Lorentz factor F^, ~ ty/(\ — Pi) ~ 2r|i^. We will take the variability scale to be = 0.01 s. 
Although variations in GRB spectra down to milliseconds have been observed, what is important here is the 
Fourier strength on the variability, and the choice can be debated. In any case, when obtaining our final 
results we will return to this problem and vary this as well as other critical benchmark parameters over a 
wide range. 

Long-duration GRB are beamed; their observed isotropically-equivalent 7-ray luminosity is on average 
Lj = 10^^ erg/s. In the fireball model it results from synchrotron radiation by electrons accelerated in 
internal shocks. The spectral photon density of the burst (GeV~^ cm~'^) in the observer's frame can be 
adequately parametrized as [28] 



,'(e/eo)"e-^/^° e < {a - f3)eo 

- /3)"-^e^-"(e/eo)^ e > (a - /3)eo ^ ' 



with a ~ —1, P ~ —2.2 and eo ~ 1 MeV. The normalization is given by = J deen^{e) — L^/Anrf. The 
corresponding photon density in the comoving frame is then given by n'^it') ~ 7i-y(F,;e'). 

After Fermi acceleration the electron population follows a power-law spectrum with minimum energy 
mill — (£e/Ce)wp in the comoving frame. Here Ce < 1 is the fraction of electrons that achieve equiparti- 
tion [5]. This spectrum yields by synchrotron radiation a photon spectrum peaking at 

(2) 

where = L^^gjlO^^ erg/s, £3 = e_B -lO.l, Eg = _i0.1, Ce = Ce -lO.l, F^ = F^^s.slO^'^ and i„ _2 = 
i^^_20.01 s. This is close to the observed peak of the GRB burst spectrum at ©(MeV). 

This concludes our description of the GRB fireball and its electromagnetic spectrum. There is noth- 
ing new here. We will next discuss the production of cosmic rays and neutrinos by applying the fireball 
phenomenology just described. 




3. Cosmic Ray and Neutrino Emission 

Although it is straightforward to argue that GRB fireballs represent an environment that can yield CRs 
of very high energy, it is important to take into account that the acceleration competes against energy loss 
due to synchrotron radiation and pion production. In other words, the fireball phenomenology is subject 
to the conditions that (i) the time to accelerate protons to the highest energy does not exceed the lifetime 
of the fireball and (ii) that the energy gained is not lost to synchrotron radiation and pion production. We 
discuss these constraints sequentially. 

For efficient acceleration of UHE CR protons their gyroradius must be contained within the acceleration 
region which is related to the size of the shock. In the comoving frame the Larmor radius of a proton is 

= E' /eB' and the corresponding acceleration time c/rj^ is given by i^^.^. = r/r^ with 77 > 1 or 

(„2^ p2 J.2 p4 \ 

where Ep — £'p,20.5l0^"'^ eV, corresponding to the upper end of the UHE CR spectrum. The size of the 
accelerator is set by the size of the shock, t'^^^ ~ ri/2Ti ~ tyVi, or 

<^y„~ 1 X 10"F,^2.5t„,-2cm. (4) 
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From this we derive the maximal proton energy in the observer's frame for which t'^y^ — iaco 

^ 4 X 10- ( X"'?^ ] ' «V (5) 

We next consider the energy losses that compete with the acceleration. In the comoving frame the time 
scale associated with synchrotron radiation is t'^y^^ = 'd'Krn^ / E' B'"^ or 

isync =i 7 X 10 "J — cm . (6) 

\ £B,-1-C'P,20.5-'j7,52 / 

The maximal proton energy is reduced when t'^y^ = t'^^^ or is smaller; in this case, using (jsj) 

i^„^3xlO- r%-^'--^f--M eV. (7) 

The photo-pion energy loss rate is determined by the pj cross section, the photon density riy and the 
average energy loss of the protons (xp^) in each interaction. In the A^-resonance approximation 

where ~ 120 MeV, = 1232 MeV and crA+ — 420 /ib and {xp^^+) ~ 0.2. Hence, for a spectral 

slope a ~ — 1 the optical depth becomes constant above a proton energy 

Ep^h ^ 0.4^GeV2 ~ 4 X lO^^^iMeV , (9) 

Co £0,6 

with eo = Co. 6 MeV. Photo-pion production therefore introduces a break in the neutron production spectrum 
at (|9]). Above this energy the energy loss length in the comoving frame is 



A inl2 / ^".6rf 2.5*1-2 \ /,„n 

4 X 10 ■ ■ — cm, (10) 



L. 



7,52 



and the maximal proton energy satisfying t'^ — t'^^^ is 



' K T \ 1/2 

22 / i,2.5'-v,-2 



Ep,^,^ 2 X 10^^ eV. (11) 



For GRB fireball parameters considered in this work the scale of photo-pion losses ( 10 1 is always larger than 
the synchrotron scale ([6| or the dynamical scale ([4|. The maximal proton energy is hence given by the 
smaller of Eqs. ([S]) and ^7]). 

There is an additional consideration as mentioned in the introduction: because protons are magnetically 
coupled to the expanding fireball they will lose energy adiabatically if they remain confined in the expanding 
shock. Adiabatic cooling, however, is not important for neutrons that are produced at scales {xp^^)t'^ and 
can escape the source |l3j . 

We will assume that the spectrum of observable CRs results from neutrons escaping the source. The 
competition between various energy loss mechanisms in internal shocks can introduce various features in the 
injection spectrum of individual GRBs |29j . Within our approximation of a homogeneous source distribution 
we assume that the effective CR emission rate can be approximated as 



{E/Ep,^)-' 



For E ^ Ep^ti, this reduces to the typical QqE^'' exp{—E / Ep „iax) approximation of the CR injection 
spectrum that we are going to test against HiRes data in the following. If energy loss of pions and muons 
prior to decay were negligible, we can relate the neutrino (per flavor) and CR neutron emission rates as 



QliE,) ~ -QcRiEje) 



(13) 



where e = (E^/En) ~ 0.06. However, synchrotron losses of secondary pions and muons in the background 
magnetic field are important. With muon and pion lifetimes of t'^'^'^ — 2.2 /is and t'^'^'^ — 26 ns, respectively, 
this will introduce a synchrotron break in the respective spectrum at 



E'-K/fJ^S 



(14) 



The corresponding break in the neutrino spectra from tt^ 



fj-^i'f_i and /i^ 



K,s ^ KlJ^ or 



E„ 



' p8 j-2 

1.2.5'^ 



1/2 



£b.-iL. 



7,52 



2xlGi7eV (z/^), 

1 X 10^6 cV (j>^,Z/e)- 



The total diffuse neutrino flux can hence be approximated as 



e+f^z/e, respectively, is 



(15) 



(16) 



where the sum runs over the neutrino flavors — or . 



4. DifTuse Cosmic Spectra 

For the calculation of the spectrum of UHE CR protons we assume that the cosmic source distribution 
is spatially homogeneous and isotropic. The comoving number density Yi — ni/{l + z)^ of nuclei of type i 
is then governed by a set of Boltzmann equations of the form: 

= dE[HEYi) + dsihY,) - r, + ^ / AE, 7,^,y, + A , (17) 

together with the Friedman-Lemaitre equations describing the cosmic expansion rate H{z) as a function of 
the redshift z: H'^{z) = [riin(l + zY + ^^a], normalized to its present value of _ffo ^ 70 kms~^ Mpc~^, 
in the usual "concordance model" dominated by a cosmological constant with JIa ^ 0.7 and a (cold) matter 
component, Q,^ ~ 0.3 [HIP. The time-dependence of the redshift is dz — —dt{l + z)H . The first and 
second term in the r.h.s. of Eq. (17) describe continuous energy losses (CEL) due to red-shift and e+e" pair 
production [31 on the cosmic photon backgrounds, respectively. The third and fourth terms describe more 
general interactions involving particle losses {i — >■ anything) with interaction rate F^ , and particle generation 
j — >■ i. We use the Monte Carlo package SOPHIA [32] to calculate P7 interaction rates and spectra. The last 
term, Ci, accounts for the CR emission rate per co-moving volume. 

We describe the cosmic evolution of the CR sources by the ansatz £cr(^, E) = 'Hgrb(2)Qcr(£')- Given 
their supernova association it is natural to assume that the comoving density of GRBs follows the star 



formation rate (SFR). We will use in the following the approximation [33l 1341 



\l + zf-^ z<l, 
Hsfr(z) - <j iVi (1 + z)-"-3 l<z<4, (18) 
\NiNi{l + z)~^-^ z>4. 
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with normalization factors, iVi = 2^-^ and iV4 = 5^'^. It has been suggested that the GRB population may 
not directly follow the SFR and may have been stronger in the past. This possibility is interesting also in 
the context of UHE CR models: a strong evolution of CR proton sources can account for the spectrum of 
extra-galactic UHE CRs even below the ankle [5S] . These "low-crossover" models predict significantly larger 
neutrino fluxes from the sources of CR protons and are already limited by upper bounds on diffuse neutrino 
fluxes [321 [37] . For the illustration of the effect of a strong GRB rate at high redshift wc will use in the 
following a comparatively strong evolution of the from [38 

Hstrong(;2) = (1 + zf'^ HsFr(^) • (19) 

This is consistent with a more recent study [39]. However, the sample of high-redshift GRBs where the 



enhanced evolution (19) becomes strongest is small and plagued by systematics. It has also been argued 
that the evolution at low redshift is consistent with the SFR whereas at high redshift it can be approximated 
as an almost constant rate [lOHH]. 

The density of CR sources at high redshift has two effects on the analysis. Firstly, a stronger evolution of 
the CR sources requires in general lower values of the power-law index 7 to reproduce the CR data. Secondly, 
a higher density of sources in the past tend to produce larger energy densities of secondary neutrinos and 
7-rays. However, we will see later on that the contribution of UHE CR proton sources to the diffuse extra- 
galactic 7-ray background is already close to maximal if we assume a source evolution following the SFR. A 



stronger evolution as in Eq. ( 19 ) can not significantly enhance the prompt neutrino fluxes without violating 
the 7-ray bound. We will estimate the effect of source evolution on the bolometric electro-magnetic energy 
density in the following. 

Particles with electromagnetic (EM) interactions produced in association with the cosmic rays, 7-rays, 
electrons and positrons, will cascade in the universal photon background and magnetic fields on time- 
scales much shorter than their production rates. The relevant processes with background photons are 
inverse Compton scattering (ICS), -I- 7bgr — > -t- 7, pair production (PP), 7 -I- 7bgr + e~ , double 

pair production (DPP) 7 + 7bgr -> e+ + + e+ + e^, and triple pair production (TPP), -I- 7bgr 

+ -\- e~ |311 . High energy electrons and positrons can also lose energy via synchrotron radiation 
on the intergalactic magnetic field with strength limited to be ~ 10~^G |44) . 

The evolution of the diffuse 7-ray and spectra follows the Boltzmann equations ( 17 ). For recent studies 



see [IHl • After cascading the EM flux accumulates into 7-rays of GeV-TeV energy with a characteristic 
and essentially universal spectrum. Its normalization is determined by the total energy density of EM 
radiation from the propagation loss of CR nuclei. For this purpose, we define the comoving energy density 
at redshift z as 

Wcas(^) = JdEE J2 Y^{z,E), (20) 

i— 7,e='= 

which fohows the evolution equation 

Wcas + ift^cas = / E)Yck{z, E) . (21) 



The energy loss rate b = dE/dt was already defined by Eq. (17), but here comprises the combined energy 
loss of CR protons into EM radiation by Bethe-Heitler production and photo-pion interactions. The energy 
density (eV cm~"^) of the electromagnetic background observed today is therefore given by 

dt fdE^Jf^ YcR(z,E). (22) 
J (1 + ^) 



Typically the energy density (20) today obtained by a detailed calculation of the EM spectra agrees with 



the (quicker) bolometric calculation ( 22 ) within a few percent and wc will use the latter in our statistical 
analysis. 

The influence of cosmic evolution on the energy density of the cascade can be estimated in the following 
way. The UHE CR interactions with background photons are rapid compared to cosmic time-scales. The 
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Figure 1: Results of the goodness of fit test of the HiRes data |48l . We show the 99% (magenta) confidence levels of the injection 
index 7 and the maximal energy i?niax- In this analysis we consider only minimal 7-ray production coming from Bethe-Heitler 
pair production and photo-pion production of CR protons during propagation ("Min"). The dashed line ("Max") assumes that 
also 7-rays from tt" production survive the internal shock environment and cascade during propagation in the inter-galactic 
background. 



energy threshold of these processes scales with redshift as £'th/(l + z) where Eth is the (effective) threshold 
today. We can therefore approximate the evolution of the energy density of the secondaries as 

Wcas + i^CJeas ^ r;casH(z) J dE E Qck{E) , (23) 
Bth/(l+z) 

where r^cas denotes the energy fraction of the CR luminosity converted to the electromagnetic cascade. 
Assuming a power-law injection Qck{E) cx E"'^ with sufficiently large cutoff i^inax S> £'th we obtain that 
strong cosmic evolution (191 enhances the diffuse 7-spectrum as Wcas oc /di(l -f z)^ '*^'''"^. For the proton 
spectrum 7 ~ 2.3 this corresponds to a relative increase of ~ 4, which agrees with numerical results. 

Photo-hadronic interactions in internal shocks produce EM radiation from the decay of pions. Most 
of this additional EM radiation is expected to cascade in the high background densities of 7-rays or 
and strong magnetic fields of the fireball and should eventually contribute to the luminosity of the burst. 
However, it is conceivable that part of it decouples from the fireball to contribute an additional source term 
Qem [l2] to the diffuse 7-ray background. Assuming energy conservation in the cascade this additional 
contribution is smaller than 



^cas. source 



Ho 



dEEQEuiE) (24) 



with ^/Hq = Jdt'H{z)/{l + z). Assuming the GRB evolution (18 1 and (19) this gives ^sfr ^ 2.4 and 
^strong — 7.3, respectively. 

In the following we will assume two limiting cases for the contribution of GRBs to the diffuse GeV-TeV 
background. The minimal model ("Min") assumes that the only contribution to the cascade results from 
the propagation from CR protons. In a maximal model ("Max") we assume that additionally all 7-rays 



from neutral pion production in the GRB contribute maximally to the cascade in the form of the term ( 24 1 
with Qem = QttO • 



5. Goodness of Fit Test 

In this section we present the details and results of our statistical analysis. We perform a goodness of 
fit (GOF) test of the compatibility of cosmic ray data with a given model, characterized by the spectral 
index 7 and the maximal energy i?max of the injected cosmic ray and by the two evolution models ( 18 ) and 
(19 1. We use CR data from HiRes I and II [H] above the "ankle" at 4 EeV. We additionally impose the 



7 



lo- 



co 

7 

<N 

I 

E 
o 

O 10 



GOF 99% E>4 EeV 

LSFR Evolution * HIRes II 
• HiRes I 
Auger 



GOF 99% E>4 EeV 

-Strong Evolution *■ HiRes II 
^ - HiRes I 

Auger 




E(GeV) 



E(GeV) 



Figure 2: The allowed proton range of proton flux at tfie 99% confidence level. The color-coding is as in Fig. [l] Each fit of the 
proton spectrum is marginalized with respect to the experimental energy uncertainty and we show the shifted predictions in 
comparison to the HiRes central values |48j . For comparison we also show the Auger data |50j which has not been included in 
the fit. 



requirement that the electromagnetic flux accompanying the cosmic rays does not exceed the Fermi-LAT 
measurements of the diffuse extra-galactic 7-ray background in the GeV-TeV energy range. 

Given the acceptance Ai (in units of area per unit time per unit solid angle) of the experiment in bin i 
centered at energy Ei with bin width and energy scale uncertainty cr^^, the number of events expected 
in the bin is 



dE J\ 



(25) 



-Ei(l-(-(5)-Ai/2 



where J^f^^^E 



[E) = np(0, E)/Att is the proton flux arriving at the detector. It is determined by the source 
luminosity of Eq. (12 1 with Ep i, ^ i?min = 4 EeV and cosmic evolution given by Eqs. (18) or (19 1. The 



parameter 5 in Eq. (25 1 is a fractional energy-scale shift that takes into account the uncertainty in the 



energy-scale and M is the normalization of the proton source luminosity. 

The probability distribution of events in the i-th bin follows a Poisson ditribution with mean Ni. Cor- 
respondingly, the r-dimensional probability distribution for a set of non-negative integer numbers k = 
{fci, ...fcr}, Pj:{n,-f,J\f,5) is just the product of the individual Poisson distributions with r the number of 
bins with Ei > E^nin- Given a model, the experimental result 7V°'^p = {A^^'^'', A'^'^p} has a probability 
Pj^cxpi'^i E^s^x,J^ , which, after marginalizing over the uncertainty in the energy scale and normalization, 
is given by 

Poxp(7: -E^max) = max^^TV' (^^--p (T: -Emax, A/", S)) . (26) 

Here the maximization is for some prior for S and Af. As a prior the energy shift S we used a top hat of 
width aE. ■ 



For Af we impose two priors. The first is associated with the upper bound on the total EM of Eq. (22 1 
that should not exceed the Fermi-LAT measurements [27], or following Ref. [45] 



(AA,7,i;n,ax) < 5.8 X 10-"cV/( 



cm" 



(27) 



The second prior on the normalization is imposed by requiring that the proton spectra do not exceed the 
HiRes I and II data below -Emin by more than three standard deviations. 
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The marginalization in Eq. ( 26 1 also determines A/bcst and (5bcst for the model, which are the values of 
the energy shift and normalization that yield the best description of the experimental CR data, subject to 
the constraint imposed by the Fermi-LAT measurement. 

In the end the model is compatible with the experimental results at a given GOF if 

^'fc(7,^^max,AAbest,<Sbost) < GOF. (28) 

Technically, this calculation is performed by generating a large number A'rcp of replica experiments following 
the probability distribution and by imposing that the fraction F with Pj: > Poxp satisfy F < GOF. 

With this method we determine the range of values of (7, i?max) that are compatible with the HiRes I 
and HiRes II data [H]. We show in Fig. [T] the regions with GOF 99% for the two evolution models ([l8| (left 



panel) and (19) (right panel). In order to illustrate the impact of the constraint imposed by the Fermi-LAT 
measurements, we also show the GOF regions without imposing it. The region bounded by the dotted line 
shows the reduced parameter space resulting from the condition that secondary EM radiation from proton 
propagation does not exceed the Fermi-LAT measurement ( "Min" ) . This bound can become even stronger 



if we consider UHE 7-ray emission from tt'^ decay in the source with diffuse energy density ( 24 ) as indicated 
by the region bounded by the dashed line ( "Max" ) . Figure [2] shows the range of the corresponding proton 
fits to the data. 

We have not included in the analysis the results from the Auger Collaboration [121 [SD] , which are shown 
in Figj2]for illustration only. As described in Refs. [351 [SO], besides the energy scale uncertainty there is 
also an (energy-dependent) energy resolution uncertainty which implies that bin-to-bin migrations influence 
the reconstruction of the flux and spectral shape. Since the form of the corresponding error matrix is not 
public, this data [49l |50j cannot be analysed outside the Auger Collaboration. 



6. Results and Their Dependence on Model Parameters 

Our final results are summarized in Fig.[3]that shows the range of prompt neutrino spectra corresponding 
to the UHE CR proton spectra at the 99% C.L. of the GOF test with the HiRes data assuming SFR evolution 
(left panels) and the stronger evolution of Eq. ( [I9| ) (right panels). We have assumed that proton acceleration 
in the shocks is efficient (77 ~ 1) and have fixed the maximum eq of the prompt 7-ray emission (Eq. (l2|) 
at 1 MeV. For each GRB model the maximum proton energy i?max corresponds to the smaller of Eq. T5| 
or Q, i.e. tdyn < tsyn Or tgyn < ^dynj respectively. For given values of fi'max, Ti and i„ we can in general 
derive two solutions of {eb /£e)L-y, which we allow to vary within the range 0.1 < {eb / £e)L^f,^2 < 10. The 
corresponding branches are indicated in the plots. 

The upper three rows of Fig. [s] shows prompt neutrino spectra with a Dopplcr factor F,; = 10^'^ of the 
internal shock. The top panels show the results for the variation time-scale = 0.01. This includes the 
range of neutrino fluxes for typical "benchmark" values of the GRB environment with = 10^^ erg/s 
and i?max = 10^°'^ eV, which is shown separately as the cross-hatched area. The range of prompt neutrino 
spectra for the benchmark model corresponds to the range of the power-law index 7 at the 99% C.L. and 
the normalization of each particular model. Note, that for both cases, SFR or strong evolution, the range 
of the prompt neutrino flux exceeds the present limits of IC-40. 

The energy density of the prompt neutrino spectrum peaks at an energy E ~ eiS'p.b — 1 PeV and is 
marginally consistent with recent upper limits on the diffuse neutrino flux by IC-40 [51) 1521 . Only for very 
short variation time-scales of ~ 10""^ s (third row of Fig. [3| and for synchrotron-loss-dominated GRBs 
(^syn < idyn) Can the predicted flux avoid present upper limits. Note that since cx r|/eo we can shift the 
peak to higher energies assuming a stronger Doppler boost of the internal shocks. However, the condition 
Ep^b ^ -E'min imposed for the fit requires that F^ < lO"^'^ and hence eEp^b < 100 PeV. 

The bottom row of Fig. [s] we show the prompt neutrino spectra for Fj = lO'^ and the full range of 
variability 10~^ < ti,/s < 0.1. Due to the strong dependence of the synchrotron time scale on the Doppler 
factor, tgy„ c>c Fj, the GRB models at the 99% C.L. correspond to proton sources where the maximal energy 
is limited by the dynamical time. The range of the flux is consistent with present IC-40 limits. Note, 
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Figure 3: Prompt neutrino spectra for SFR (left) and strong (right) evolution corresponding to the 99% C.L. of the "Min" 
model shown in Fig. [l] and assuming the luminosity range 0.1 < (£s/ee)i'7,52 < 10. We show the prompt spectra separately 
for the branches fjyn < isyn (green right-hatched) and tsyn < tdyn (blue left-hatched). The preliminary IceCube hmits |51II52I 
(90% C.L.) on the total neutrino flux from the analysis of high-energy (HE) and extremely-high energy (EHE) neutrinos with 
the 40 string sub-array (IC-40) assume flavor composition after osciflation. We also show the sensitivity of the full IceCube 
detector (IC-86) to EHE neutrinos after three years of observation. The gray solid area shows the range of GZK neutrinos 
expected at the 99% C.L. 
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however, that the opacity of of the GRBs for UHE p7 interactions and subsequent emission of CR neutrons 
is very sensitive to the Doppler factor, 



Increasing the average Doppler factor of the GRB from 10^'^ to 10"^ reduces the efficiency of UHE CR 
emission by two orders of magnitude. It is hence more hkely that the successful candidates of CR sources 
are GRBs with lower F^. 

The range of prompt neutrino spectra shown in Fig. [3] spans about two orders of magnitude, a range 
larger than a previous study in Ref. 53 . The difference of the results can be traced back to different 
assumptions. Firstly, we here assumed that only neutrons can escape the GRB shock environment and form 
the spectrum of UHE CRs. Hence, both, CR and neutrino spectra depend on the opacity of the source. 
Secondly, we apply a statistical fit to the data assuming a wide range of spectral indices. Statistically allowed 
CR spectra with a steep injection index 7^2 correspond to GRB models with a higher bolometric energy 
density and hence larger prompt neutrino emission. For large maximal proton energies, -Emax S> 10^^ eV, 
we can approximate the total local power density in UHE CRs derived from the best fit to the HiRes data 
as 

ecR= / dSS/:cR(0,i;)^(l-2)xl044ergMpc-3yr-ix <^ V Y ^ ' (30) 

This value is consistent with the original study by Waxman and Bahcall jl5) assuming a total power density in 
UHE CRs of the order of 4 x lO^* erg Mpc"^ yr'^ between lO^^ eV and lO^i eV assuming a fiat spectrum with 
power-law index 7 = 2. In our approach we allow steeper injection spectra 7 3> 2 extending down to the break 
energy Ep^\, <^ 10^^ eV (cf. Eq. ([o])) which is set by the GRB fireball environment. For very steep injection 
spectra with 7 ~ 2.6 consistent with our fit at the 99% C.L. (cf. left panel in Fig. [l| and i?p_b = 4 x 10^^ eV 
this corresponds to a maximal local power density of the order of ecR.max — lO'^^rg Mpc~^ yr~^. 

These large local power densities of steep UHE CR proton spectra are also consistent with the results 
of Ref. contrast to Ref. [51] we include here the limits implied by the extra-galactic diffuse 7-ray 

background inferred by Fermi LAT. This bound has the strongest effect on strong evolution models with 
large power-law index as can be seen in the right panel of Fig. [l] The corresponding proton fluxes at the 
99% C.L. are shown in the right panels of Fig. [2] and the prompt neutrino fluxes in the right panels of 
Fig.|3] The range of spectra is comparable to the moderate evolution following SFR (left panels) since the 
contribution of UHE CRs to the diffuse extra-galactic 7-ray background is already close to maximal. 



7. Conclusion 

We have discussed prompt neutrino emission associated with the production of UHE CR protons in 
GRBs. Our analysis assumes that UHE CRs above the "ankle" at 4 EeV consists of neutrons emitted from 
P7 interactions in internal shocks of the GRB fireball model. We determined the CR emission density by a 
fit to HiRes data and used this information to determine the corresponding neutrino fiuxes for a wide range 
of parameters for the shock environment. We have also carefully studied the effect of secondary production 
during CR propagation in the form of GZK neutrinos as well as the GeV-TeV 7-ray background. 

The main results are shown in Fig. |3] We have shown that typical ( "benchmark" ) fireball environments 
predict prompt neutrino fiuxes that exceed present diffuse neutrino limits from IceCube (IC-40) if the 
associated proton spectrum is fitted to actual CR data. This remains partially true if we relax the condition 
of internal shock parameters. Consistency with the diffuse neutrino limits requires that the CR acceleration 
takes place in GRB fireballs with small internal shock radii (corresponding to variation time-scales of milli- 
seconds) and/or large Lorentz boost ~ 1000. Alternatively, it is possible that UHE CRs at the ankle 
might still receive a significant contribution from galactic CRs, which lowers the required power density in 
extra-galactic CRs and the corresponding prompt neutrino spectra from GRB sources |55j . However, one 
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has to worry that galactic CRs at these energies will reveal themselves via an anisotropy toward the galactic 
center. In any case, the sensitivity of the full IceCube observatory (IC-86) after three years of observation 
will improve present diffuse neutrino limits by a factor five and will serve as a crucial test of the GRB 
scenario of UHE CRs. 

Finally, we would also like to mention additional constraints of this CR scenario. A fraction of UHE 
7-rays produced in neutral pion production may escape the fireball shocks and contribute to the GeV-TeV 
background as well. As an estimate, we have tested a pessimistic scenario ( "Max" model) where all 7-rays 
from 7r° production escape the source. We found that this additional contribution may become important 
in the case of source evolution much stronger than the SFR (c/. dashed regions in Figs. [T]and[2]). It is 
also possible to constrain UHE CR emission from internal shocks of GRBs by the integrated fluence of all 
bursts as has been done recently in [56 . On the other hand, the contribution of the burst spectrum to the 
diffuse 7-ray background from unidentified cosmic GRBs has been discussed in [57) and has been shown to 
contribute a relatively small fraction. 
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